Absence of a Direct Superfluid to Mott Insulator Transition in Disordered Bose 

Systems 
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We prove the absence of a direct quantum phase transition between a superfluid and a Mott 
insulator in a bosonic system with generic, bounded disorder. We also prove compressibility of 
the system on the superfluid-insulator critical line and in its neighborhood. These conclusions 
follow from a general theorem of inclusions which states that for any transition in a disordered 
system one can always find rare regions of the competing phase on either side of the transition 
line. Quantum Monte Carlo simulations for the disordered Bose-Hubbard model show an even 
stronger result, important for the nature of the Mott insulator to Bose glass phase transition: The 
critical disorder bound, A c , corresponding to the onset of disorder-induced superfluidity, satisfies 



the relation A c 



> -Eg/2, with -Eg/2 the half- width of the Mott gap in the pure system. 
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PACS numbers: 03.75.Hh, 67.85.-d, 64.70.Tg, 05.30.Jp 

The interplay between disorder and interactions is one 
of the most exciting and longstanding problems in con- 
densed matter physics, full of controversial theoretical 
and experimental results. Bosonic systems are especially 
difficult to handle theoretically because the limit of van- 
ishing interactions is pathological: in the ground state all 
particles occupy the same lowest-energy localized state. 
Hence, interactions and disorder are equally important 
at the quantum phase transition between superfluid (SF) 
and insulating (I) ground states. Relevant examples from 
Nature include 4 He in porous media and aerogels 111, 4 He 
on substrates 0, Q , thin superconducting films (3[ , and 
Josephson junction arrays [1]. 

The physics is further complicated when bosons are 
subject to a periodic external potential and interactions 
are strong enough to drive a commensurate system (with 
particle number an integer multiple of the number of lat- 
tice sites) to the insulating Mott phase (MI) even without 
disorder. A fundamental problem here is the role of an 
arbitrarily weak disorder in the vicinity of SF-MI quan- 
tum critical point of the pure system. 

Building on one-dimensional results by Giamarchi and 
Schulz [j| Fisher et al. argued the existence of a Bose 
glass phase (BG) in the disordered Bose-Hubbard Hamil- 
tonian in any dimension [|| : while the commensurate MI 
is gapped, the insulating BG remains compressible. Ever 
since, the Bose-Hubbard model has received a lot of the- 
oretical attention, although it remained beyond direct 
experimental reach. This is changing with the experi- 
mental demonstration of the SF-MI transition of ultra- 
cold atoms in optical lattices 0] following the theoretical 
proposal of Ref. Q. At present, cold- atom systems of- 
fer an unprecedented control over system properties and 
have been explicitly shown to be accurately described 
by the Bose-Hubbard model. In particular, there is full 
agreement between the experimental data and Quantum 
Monte Carlo simulations, proving that the optical lat- 
tice system can be a quantitatively reliable simulator Q . 



Experiments have recently advanced to the stage where 
disorder can be added and controlled using speckle po- 
tentials [IS III III. 

Though the pure system is well understood, the 
properties of the disordered Bose-Hubbard model 
remain subject to debate, even at the qualitative 
level. The relevance arbitrarily weak disorder at the 
superfluid-insulator quantum phase transition at an 
integer filling factor is a highly controversial issue 
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19|, [20J, |21|, 122, 

Fisher et al. 



raised the 



problem and argued that a direct SF-MI transition was 
unlikely (see also Refs. [3ll |32|). though not fundamen- 
tally impossible. Curiously, a large number of direct 
numerical simulations fl5l. Il8l. [TpL pll. 26l. 27l and some 
approximate approaches observed 
this unlikely scenario! 

However, all numerical simulations reporting a direct 
SF-MI transition were ignoring a rigorous theorem [|| [H[ 
(to be referred to as Theorem I) stating that if the bound 
A on the disorder strength is larger than the half-width 
of the energy gap E g / 2 in the ideal Mott insulator, then 
the system is inevitably compressible, i.e. the transition 
is to the BG insulator and not the MI whenever 



A c > E. 



/2- 



(1) 



Intuitively, the condition A > E g / 2 (U) seems to be 
necessary for the disorder to at least destroy the Mott 
gap, putting aside the question of onset of superfluidity 
which is likely to require even stronger disorder. There- 
fore, it is reasonable to conjecture that Eq. (pj holds for 
the Bose-Hubbard model (see Refs. [HI, [32[ for more de- 
tails). This condition was already shown to be necessary 
in ID 22], and there is compelling numerical evidence 
that it equally holds in 2D [28|. Nevertheless, a rigorous 
proof valid for any dimension has not been found. More- 
over, it was recently claimed [29| that the relevance of 
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weak disorder depends on dimensionality and that in 3D 
weak disorder can be irrelevant. 

In this Letter, we introduce and prove the theorem of 
inclusions: In the presence of generic bounded disorder 
there exist rare, but arbitrarily large, regions of the com- 
peting phase across the transition line. By generic disor- 
der we mean that any particular realization has a non- 
zero probability density to occur in a finite volume. This 
theorem immediately implies the absence of a direct SF- 
MI quantum phase transition. 

We will also introduce Theorem II which states a non- 
zero compressibility on the superfluid-insulator critical 
line and in its neighborhood for models which have dis- 
ordered on-site potentials. Theorem II is not based on, 
and thus does not imply condition ([T]) . To check whether 
Eq. fl]) holds in three dimensions, as well as to produce 
some quantitative benchmarks for future experiments, we 
complement the theorem by first-principles simulations 
of the disordered Bose-Hubbard model by the worm al- 
gorithm [lH . For weak disorder, our results are generic 
through standard long-wavelength universality consider- 
ations. 

The Hamiltonian of the Bosonic Hubbard model with 
on-site disorder reads: 

The subscripts label the sites of a (simple cubic) lattice, 
a\ is the boson creation operator, the symbol de- 
notes nearest-neighboring sites, rii = a\a.i is the density 
operator. In what follows we use the amplitude of the 
hopping term t as the unit of energy, U is the on-site pair 
interaction, fi is the chemical potential, and £j represents 
the random disorder potential. Normally, one assumes — 
without loss of generality — that there are no correlations 
between the ej's on different sites, and that their values 
are uniformly distributed on the [—A, A] interval. 

The proof of Theorem I is based on the fact that in the 
infinite system one can always find arbitrary large 'Lif- 
shitz' regions where the chemical potential is nearly ho- 
mogeneously shifted downwards or upwards by A. There 
is no energy gap for particle transfer between such re- 
gions, and they can be doped with particles/holes [34| . 
Note also that Theorem I immediately implies that the 
MI phase does not exist in a system with unbounded 
disorder. 

This proof gives a hint of a potential pitfall when one 
tries to directly interpret numerical or experimental data 
on finite system sizes. For A — > the distance be- 
tween regions contributing to a non-zero, but exponen- 
tially small, compressibility is diverging exponentially. 
If disorder is not very strong, a typical finite-size clus- 
ter is not supposed to contain even a single rare region, 
and the system behavior perfectly mimics a direct SF-MI 
transition. At the same time, Theorem I also provides a 
straightforward way out. One just needs to compare the 
critical amplitude of the disorder, A C (U), to E g / 2 (U): If 



condition (JTJ) is satisfied, then the transition is towards 
the Bose glass rather than the Mott insulator, with the 
compressibility remaining finite at the critical point. We 
now proceed to the proof of the theorem of inclusions. 

Proof of the theorem of inclusions. The proof is 
straightforward if reducing the disorder bound enhances 
the supcrfluid phase (cf. e.g. Refs. [H, [HI)- In the op- 
posite case, let us introduce £ denoting all microscopic 
parameters other than the bound A that fix the shape 
and correlations of the disorder distribution. In the pa- 
rameter space (U, A, £ ), the critical hyper-surface of the 
supcrfluid to insulator quantum phase transition is writ- 
ten as A = A C (U,£). For simplicity, we keep U fixed, 
without loss of generality. On physical grounds, A c (£ ) is 
a continuous function of its microscopic parameters, and 
a generic £ is not an extremum of A c . We do not require 
that A c (£) is analytic, even though this is most likely 
to be the case physically. Generic disorder implies that 
there exist arbitrarily large regions in which the disorder 
realization generated by £ can be considered with finite 
probability density as a typical realization of a different 
set £fl U ct- In a sufficiently small neighborhood of the crit- 
ical surface, we can always find £fl uc t such that all points 
A G [A c (£ ), A c (£fl uc t)] lie in a supcrfluid domain of the 
phase diagram thanks to £ being non-extremal. The de- 
viation from the critical surface, SA = A c (£fl uct ) — A c (£ ), 
can be chosen small enough to guarantee that it can be 
compensated by changing £ to £fl uc t- We conclude that 
within a distance \8A\ from the critical surface the sys- 
tem develops arbitrarily large supcrfluid domains, which, 
by virtue of standard Lifshitz tail arguments implies the 
absence of the gap (as in Theorem I) and thus the absence 
of a direct SF-MI transition. 

Proof of theorem II. Close enough to the critical sur- 
face we can choose a small finite value 5A < since 
fluctuations of £ can always be used to compensate for a 
sufficiently small deviation from the superfluid domain. 
In such a rare superfluid domain, the chemical poten- 
tial can be homogeneously shifted by a finite amount 
Sfi ~ 5A without violating the global disorder bound. 
Since one can always find rare supcrfluid regions large 
enough such that the finite-size quantization of the en- 
ergy is smaller than the allowed variation Sfi, the density 
of states at the global chemical potential is guaranteed to 
be finite, i.e. the system is compressible across the SF-I 
transition. Note that our proof is not based on a specific 
form of the Hamiltonian and thus applies to all bosonic 
systems with the on-site potential disorder featuring SF-I 
quantum phase transition. 

Numerical analysis. In our numerical study of model 
([2]) we consider Si uniformly distributed on the [—A, A] 
interval, and choose the chemic al p otential in the mid- 
dle of the n = 1 Mott lobe [37| to ensure that the 
system is always in the (n) = 1 state. We employ 
the worm algorithm, implemented in two independent 
ways [HI, HH ■ Our finite-size scaling needed for the SF-I 
transition is based on the square of the winding number 



30.5 30.6 30.7 30.8 30.9 31 
U/t 

FIG. 1: (Color online). Winding number squared as a func- 
tion of interaction strength (U/t) for A = 5t using a scaling 
analysis with the temporal exponent z — 2, chosen for nu- 
merical convenience (see text). The corresponding pairs of 
the linear system size and inverse temperature are : L = 4 
and pt = 1, L = 8 and fit = 4, L = 16 and fit = 16, L = 24 
and pt = 36, and finally L = 32 and pt = 64. The values of 
the chemical potentials were taken from Ref. [3?]]. 
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FIG. 2: (Color online). Finite size scaling of critical points 
(U/t) c {L) as a function of 1/L 2 for disorder strength A = 5t. 
The data points shown here are the intersection points in 
Fig. [T] In the thermodynamic limit we find (U /t) c = 30.57 ± 
0.02 by extrapolating according to (1/L) 2 (see text). 



of world lines (W 2 ), related to the superfluid stiffness by 
A s = (W 2 )/TLd [Hg. In Fig. Q] we show (W 2 ) for differ- 
ent values of U jt and system sizes at fixed disorder bound 
A = 5t where we performed scaling with the temporal ex- 
ponent z = 2 which is our choice to reduce the numerical 
effort. The intersection points of the curves with dif- 
ferent system sizes, are shown in Fig. [2j Extrapolation 
to the limit of infinite system size using a (1/L) 2 scal- 
ing yields a critical value (U/t) c = 30.57(2). Although 
the exponents for corrections to scaling are not known 
precisely, the critical values (U /t) c are determined very 
accurately; we have repeated the simulations with z = 1 
(not shown) and found the same value of (U/t) c within 
error bars, showing that the choice of z does not mat- 
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FIG. 3: (Color online). Phase diagram of the 3D Bose- 
Hubbard model with disorder. The E g / 2 (U) curve is taken 
from Ref. [33] and is shown here to demonstrate that the con- 
dition A c (U) > E g / 2 (U) is definitely satisfied. It also marks 
the BG-MI transition boundary according to the conjecture. 
Error bars are shown, but are smaller than point sizes. 



ter for the location of the critical point. Indeed, winding 
numbers are exponentially suppressed in the insulating 
phase, and thus any temporal exponent is supposed to 
yield the same critical point in the scaling limit. It is 
expected that close to the tip of the Mott lobe the in- 
fluence of the U(l) symmetry is strong enough to pre- 
vent one from directly (e.g. by using critical exponents) 
distinguishing between different universality classes from 
finite system simulations. For other disorder bounds we 
find (U/t) c = 29.44(2) for A = 1.25i, (U/t) c = 29.64(2) 
for A = 2M, and (U/t) c = 30.57(2) for A = lOt. In the 
latter case, we directly observed that at the transition 
point the particle number fluctuations diverge with the 
system size. Still, even for this relatively strong disorder, 
finite compressibility can be revealed only by going to 
very large samples, L > 32, meaning that the compress- 
ibility of the Bose glass phase remains tiny and quickly 
diminishes away from the critical point. 

From our numerical data we extract the phase diagram 
of the commensurate system shown in Fig. [3] It is clear 
that the critical disorder strength is always larger than 
half the energy gap of the homogeneous model by a large 
margin. Theorem I then immediately leads to the conclu- 
sion that the SF-I transition is always to the compressible 
Bose glass phase, in compliance with Theorem II. 

In summary, we have presented an analytic proof of the 
absence of a direct SF-MI quantum phase transition in 
a bosonic system with generic, bounded disorder which 
resolves a two-decade long controversy. We have shown 
that any bosonic system with on-site potential disorder is 
compressible at the superfluid-insulator critical line and 
in its neighborhood. The quantitative analysis of the 
Bose-Hubbard model in 3D, based on worm algorithm 
Monte Carlo simulations, reveals additional details and 
provides important numerical benchmarks for future ex- 
periments. In the limit of small disorder amplitude, the 
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supcrfluid to insulator transition takes place at the in- 
teraction strength U > Uc , with being the Mott 
transition point of the pure system. The overall picture 
is that disorder first destroys the Mott insulator, con- 
verting it into a compressible Bose glass, which is guar- 
anteed to happen (and presumably happens exactly at 
the condition) when the bound on the disorder strength 
reaches the value A = E g / 2 (U), with E & /2{U) the half- 
width of the Mott gap in the pure system. The criti- 
cal bound on the disorder, A c (?7), corresponding to the 
onset of disorder-stimulated superfluidity, satisfies rela- 
tion |T]), which guarantees the above-mentioned picture. 
Thus, disorder is relevant to the SF-I criticality in all 
experimentally accessible dimensions. 

The theorem of inclusions is readily generalized to 
an arbitrary transition in a disordered system between 
phases A and B, and states that one can always find arbi- 



trarily large inclusions of the competing phase on cither 
side of the transition line. The SF-I case with diago- 
nal disorder considered here is one example; the generic 
superfluid-Mott(or checkerboard)-glass transition in a 
particle-hole symmetric system with off-diagonal disor- 
der is another one (cf. [HI, Whether rare regions 
are relevant to the nature of the phases A or B remains 
model specific and is not part of the Theorem. 
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